Transient features of quantum open maps 
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We study families of open chaotic maps that classically share the same asymptotic properties 
- forward and backwards trapped sets, repeller dimensions, escape rate - but differ in their short 
time behavior. When these maps are quantized we find that the fine details of the distribution of 
resonances and the corresponding eigenfunctions are sensitive to the initial shape and size of the 
openings. We study phase space localization of the resonances with respect to the repeller and find 
strong delocalization effects when the area of the openings is smaller than h. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

The quantum treatment of scattering situations where 
the trapped set has a complex fractal structure has re- 
cently been the source of much attention, particularly 
with respect to the confirmation of the fractal Weyl law 
[l|-[3[ hypothesis concerning the global distribution of 
long lived resonances and their localization properties on 
the classical repeller [4-6] . The use of quantum maps to 
model these systems has been of paramount importance, 
allowing extensive numerical, and sometimes, analytical 
results. 

A classical "closed" map B is represented as a diffeo- 
morphism on a symplectic manifold S [TJ. Its open ver- 
sion B can be defined restricting the phase space to V 
(a subset of S), which means that trajectories outside 
of these region escape. The associated classical repeller 
/C is defined as the set of trajectories which remain on 
V both in the past and the future for infinite times (i.e. 
the intersection of the backwards and forward trapped 
set). If the closed map B is chaotic, the dynamics of B 
on the repeller is hyperbolic and complex, and the result- 
ing /C has a fractal structure with non-integer dimension 
d. Also, an initially uniform probability distribution in 
phase space decays exponentially with the classical es- 
cape rate 7. These are asymptotic features of the map 
appearing as T ^ 00, that are associated to the global 
distribution of resonances and how their number scales 
with h. Note that given an open map we have a unique 
associated repeller, but the converse is not true. A given 
repeller can be related to a large set of open maps (e.g. 
for the same S, the repeller given by V and V" = BiV) 
are the same). 

When quantized, the evolution of these systems is non- 
unitary, there is a probability loss reflecting the clas- 
sical loss of trajectories. The corresponding operators 
have right and left decaying nonorthogonal eigenfunc- 
tions with complex eigenvalues \j , also referred to as res- 



* Email address: 'ermann@tandar.cnea. gov. ar| 



onances. They fall inside the unit circle in the complex 
plane (|Ajp = exp(— Fj) < 1). The exponent Vj > can 
be interpreted as the rate at which probability is lost and 
thus is called the (quantum) decay rate. The long lived 
resonances are those with the lowest decay rates and are 
directly associated to the classical repeller /C. Their Vj 
are related to the classical escape rate 7 in the generic 
case (see however Q). 



In this paper we study classical maps having the same 
repeller and escape rate but different openings, and their 
corresponding quantum counterparts, by means of the 
paradigmatic tri-baker map. How these differences affect 
the behavior of the resonances is a fundamental question. 
Moreover, this is a crucial step in order to understand 
and evaluate the performance of the semiclassical theory 
attempting to describe open maps ^9|] . Out of the many 
ways to open a map which lead to the same repeller we 
have found that two of them are the most important for 
this purpose, giving rise to two families of maps. Keeping 
the area of the opening constant {shift family) or letting 
it to increase {intersection family) defines them and is 
the source of important dynamical differences during the 
first time steps (though the asymptotic decay rate is the 
same for both of them). In the first case, we have found 
that by just varying the shape of the opening, the be- 
havior of the long lived portion of the spectrum changes 
slightly, in a non trivial fashion. In the second case, the 
quantum spectrum changes as a whole. Finally, the shift 
family allows to identify a quantum phase transition in 
the localization of the resonances, for openings having 
areas smaller than h. 



In the following we describe the organization of this pa- 
per. In Section [ni we define the family of maps that we 
study and their quantizations. In Section Hill we analyze 
the properties of the spectrum and the corresponding set 
of eigenfunctions by means of different measures and a 
phase space distribution that is specially useful for open 
maps [10] . Finally, in Section IIVI we point out our con- 
cluding remarks. 



II. OPEN MAPS 



[0, 1) X [0, 1) is given by 



Given a repeller /C there are many open maps that have 
it as the invariant set. We will take Vi, an opening in 
phase space, as the generator of all these open maps with 
the same repeller. For simplicity, in our example (which 
follows) it will be taken symmetric in q and p (invariant 
under q ^ p exchange). Other choices that do not re- 
spect this symmetry (or equivalent ones in other cases) 
lead to more differences in the resonances whose study 
is out of the scope of the present paper [ll| . The next 
step in the construction consists of performing the evo- 
lution of Vi under S^ that wih be called Vt = B^-\Vi). 
Finally, simple combinations of Vt for positive and nega- 
tive times allow to create the open maps, which can be 
grouped into map families according to their properties, 
all of them sharing the same repeller /C. Moreover, iter- 
ation of open maps B* can be used alternatively to build 
new maps corresponding to this invariant set. 

We will focus into two families of maps that are es- 
sentially given by preserving the area of the opening for 
all of its members in one hand or not preserving it on 
the other. These families are the most interesting since 
they allow to understand the dependence of the quantum 
spectrum on the shape and on the size of the opening, 
respectively while keeping the repeller and asymptotic 
decay rate invariant. The members of the first family 
S| with k = 1,2,..., which we call shift family, are de- 
fined as one iteration of the closed map B followed by 
the escape of trajectories outside the region 14 . In turn, 
the intersection family is defined as B^ with /c = 1, 2, . . ., 
where the allowed region this time is given by the inter- 
section of regions until time /c, Vk H Vk-i H . . . H F2 H Fi. 
For this family, the allowed region decreases for larger k 
members. Of course, for nontrivial B and Vi all family 
members are different, though the first member of both 
families are exactly the same by definition (Sf = Bl). 



The tri-baker map 

We have chosen the paradigmatic tri-baker map as the 
model for our studies. This is one of the simplest chaotic 
maps which can be easily described by a ternary Bernoulli 
shift, and where the openings can be done following sta- 
ble and unstable manifold directions. The general open 
map is defined as the composition of the closed tri-baker 
transformation followed by a given opening. We will 
define both shift and intersection families generated by 
an initial allowed region Vi having two horizontal strips 
p G [0, 1/3) U (2/3, 1). The openings wih be defined to 
be symmetrical in q and p, in such a way as to maintain 
the time reversal symmetry of the closed map as shown 
in the left panel of Fig(T] {k = 1). 

The tri-baker map in a unit square phase space T^ = 



B 



3q - [3q] 
(P+M)/3 



(1) 



where [q] means the integer part of q. The map 
is uniformly hyperbolic with Lyapunov exponent A = 
log 3. The symbolic notation of the map action is given 
by a Bernoulli shift in ternary representation of q = 
0.606162 . . . and p = 0.6_i6_26_3 . . . (given by the cor- 
responding trits 6i = 0, 1, 2) as 

(p|q) = ...6_26_i.6o6i... -> (p'lqO = ...6_26_i6o.6i..., 

(2) 
where the dot is moved one position to the right. Differ- 
ent openings of both shift and intersection families can be 
straightforwardly defined in ternary notation using open 
trits 6 with forbidden value 1 (e = 0, 2). 

Shift family members Bf^ have two open trits corre- 
sponding to the k-th most significant trit of both posi- 
tion and momentum in ternary representation. On the 
other hand, intersection family members have the first 
k open trits of both position and momentum. In order 
to illustrate shift and intersection family members, S3 
and S3 are shown in the next equation with open trits 
highlighted inside boxes 



S| > ...6-56-4 6-3 6-26-1.6061 62 6364... (3) 



Bi 



. . . 6-56-4 



e-3 6-26-1.6061 62 



6364 ... (4) 



The first 3 members of both families are represented ge- 
ometrically in Figdl In this way the classical escape 
rate can be computed analytically. When a closed trit 
is opened (transformation from Cj -^ ij ) the area of the 
allowed space is reduced by a factor 2/3 {A -^ 2/3^4). 

For the case of the k member of the shift family, two 
trits are opened in the first iterations of the map until the 
first open trit of position {ek-i) reaches the k-th most 
significant trit of momentum. From then on the map 
only opens one trit in each iteration. Therefore, during 
the first 2k iterations the allowed area decreases by a 
factor of (2/3)^ in each step, and from then the area 
decreases by a factor of 2/3. The same reasoning applies 
to the k member of the intersection family, where 2k trits 
are opened in the first iteration, and from then only one 
closed trit is opened at each step. Therefore the allowed 
area in the first step decreases suddenly by a factor of 
(2/3)^^, but from the second step decreases by (2/3) and 
so on. Thus the asymptotic escape rate is the same for 
both families. 



Quantum version 

There are different ways to quantize the tri-baker map 
[l2| . We will follow the Balazs-Voros-Saraceno (BVS) 
quantization, which is done by choosing anti-periodic 
boundary conditions in q and p [13]. The Hilbert space 
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Figure 1: (Color online) Phase space representation of the 
openings for shift and intersection families with /c < 3, in 
top and bottom panels respectively. The allowed regions are 
shown in white while the corresponding openings are in light 
blue (gray). 



dimension is given by the integer D = l/(27rfi) with po- 
sition and momentum eigenvectors \qj) and \pj) (where 
j = 1, . . . ,D) connected by the anti-symmetric Fourier 
transform Gd^ given by 



(Gd) 



1 



'\Pj) = —r^^ 



:^(i'+|)(i+|) 



D 



(5) 



The quantum baker map on a Hilbert space of dimension 
D = 3M can be achieved converting the most significant 
qutrit of position in the most significant qutrit of momen- 
tum. The matrix of the map in mixed representation has 
the form of three diagonal blocks with a finite Fourier 
transform of size D/3 in each one of them. In position 
representation we have 



B„ 



GjjBjj 




(6) 



The map is opened by means of a projector H, and in 
order to respect the q^ p symmetry we will use quantum 
open maps of the form B^i^ — H^^^a^n 

The quantum version of the shift and intersection fam- 
ilies can be easily implemented for the quantum tri-baker 
map on / qutrits with D = SK The forbidden 1 in 
one trit e, can be quantized by the one qutrit projec- 
tor tt =/ - |1)(1| = |0)(0| + |2)(2|. In this way, the 
projector applied to the i-th qutrit can be written as 



n, = / ( 



Ott I 



(7) 



with n^ = nj, nf = n^ and [n^,n^] = O where 
i, j = !,...,/. Following this notation, the shift and in- 
tersection families of open quantum baker maps can be 
straightforwardly defined as 



Bt = Gn^kB 



Bl = 



gM 



nllt 



nl 



-l^kB-mix^kJi 



fci-lfe-1 



...Hi 



(8) 
(9) 
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Figure 2: (Color online) Spectrum in the complex plane for 
three different BVS quantum open baker maps with dimen- 
sion D = 3^ = 243. Black circles represent the eigenvalues for 
the first member of both the shift and intersection families 
(/c = 1), red squares stand for shift family with /c = 5, and 
blue diamonds for intersection with /c = 5. 



with /c = 1, . . . , /, respectively. 

A phase space representation of the openings for shift 
and intersection families with /c < 3 is shown in Fig. [H 



III. RESULTS 

A. Distribution of resonances 

The quantum open maps B defined in Sec. [Ill are repre- 
sented by non- normal matrices of dimension D ([^, B'^\ ^ 
0). Therefore, their spectrum and eigenstates are given 

by 



B\Rj) 
{Lj\B 






(10) 

(11) 



where the complex eigenvalues Xj are inside unit circle 
|Aj| < 1 and \Rj) and \Lj) are the right and left eigen- 
states with j = 1,...,D. For a given eigenvalue, the 
corresponding left and right eigenstates are not equal in 
general, but they obey the orthogonality rule {Lj\Rj>) ex 
dj^j'. In Fig. [2] we show the spectrum corresponding to 
both families. It can be seen that for the first member 
of each family the results are the same as we expected 
by definition. But for the last members allowed by the 
dimension D = 3^ = 243, i.e. /c = 5 we find great differ- 
ences. In fact, the one belonging to the shift family has 
a spectrum similar to the previous case while the one of 
the intersection family is highly contractive (it shrinks 
towards zero) due to the larger area of the openings. 

In order to have a more detailed idea about the shape of 
the spectrum of both families, in Fig. [3] we show the his- 




Figure 3: (Color online) Distribution of resonances for shift 
(top panel) and intersection (bottom panel) quantum open 
baker maps with dimension D = 3^ = 2187. Moduli of eigen- 
values are represented by black circles for k = 1, 3, 5, 7 from 
top to bottom, respectively. Histograms of these moduli for 
100 intervals from to 1 are shown by means of red lines. 



tograms corresponding to the distributions of the moduli 
of their eigenvalues for I) = 3^ = 2187. 

It can be seen that the shift family has a non monotonic 
behavior with respect to the significance of the open trit. 
The long lived sector of the resonances, clearly visible as 
a big bell towards the rightmost part of the histograms, 
"oscillates" from the right to left and back, while more 
or less keeping its shape. This oscillation shows that, 
though the classical escape rate is the same for all the 
maps, the average value of the decay rate of the long 
lived states changes. 

In the case of the intersection family there is a com- 
pletely different behavior. At very short times (in fact, 
from the second time step on) the classical escape rate 
settles down to its asymptotic value. But the allowed 
area decreases suddenly also. This fact translates into 
the sudden shrinking of the spectrum as k increases. In- 
deed, the last member of this family allowed by the di- 
mension D has all its eigenvalues heavily concentrated 
around modulus zero. Then, it seems that one has to 
be aware of short time dynamics at the time of count- 



ing resonances; i.e., not only the escape rate and fractal 
dimension counts. 



B. Eigenstates and the quantum repeller 

We investigate the morphology of the eigenstates of the 
open tri-baker map families by making use of a recently 
introduced phase space representation suitable for open 
systems [10]. The first step is to define the symmetrical 
operators of right and left eigenstates 



hj 






(12) 



which are related to the eigenvalue \j and are inde- 
pendent from eigenstate normalization. These operators 
have the advantage that they clearly show the repeller 
structure underlying in the eigenstates. They are equiv- 
alent to projectors on the eigenstates of unitary evolution 
operators. To have a better idea of how the set of long 
lived resonances distribute on the classical repeller we 
also study the generalized operators onto this set taken 
as the sum of the first j of the hj' projectors, ordered 
by decreasing modulus of the corresponding eigenvalues 
(|A,|^|A,v|withj</) 



Qj = Yl ^^' 



(13) 



i'=i 



From the definition of hj in Eq. [T2]we have that tr{Qj) = 
j and Q'^j = Qj, but Q] ^ Qj. 

The phase space representation of both h and Q oper- 
ators can be defined by means of coherent states \q,p) in 
the same way as in ^10|] 



hj{q,p) = \{q,p\hj\q,p)\ 

QjiQ^p) = \{q^p\Qj\q^p)\ 



(14) 
(15) 



Clearly, the hj and Qj distributions tend to be localized 
on a discrete time version of the repeller (for Ehrenfest 
time /). Individual hj are distributed in the discrete time 
repeller in the same way as eigenstates of closed systems 
spread on the torus. Worth mentioning, scarring phe- 
nomena is visible for some hj on periodic orbits that be- 
long to the repeller (see panels in Figs. |4]and[5]). 

In order to obtain a better characterization of the 
eigenstates we quantify the localization in phase space 
by defining the norm ratio {nr) of real and positive dis- 
tributions hj (and Qj): 



nr{h) 



\\Pc\U\\Pch 



(16) 



where a coherent state is used for normalization pc 
\q^p){q^p\^ cind the phase space norm is defined as 



\h\\ 



T2 



h{q,pydqdp 



(17) 



K 




Figure 4: (color online) Norm ratio (nr) for hj operators of 
the shift quantum open baker maps with dimension D = 3^ — 
243. Black circles, red squares and blue diamonds represent 
k — l^k — 2 and /c = 5 members respectively, nr is also shown 
for corresponding Qj in dashed lines. Three Husimi functions 
(hj{q,p) = \{q,p\hj\q,p)\) for /c = 5 and j = 1, j = 13 and 
j = 28 are shown going form zero (blue) to maximum (red) 
on top panels (from left to right). Horizontal dot-dashed blue 
line is the norm ratio value obtained for the uniform quantum 
distribution on the repeller Irep- 



The norm ratio is independent of the h normalization, 
and the position in phase space (q^p) of the coherent 
state in pc. Note that norm ratio is 1 for a coherent state 
by definition, and reaches the maximum value D/2 for 
a uniform distribution in a D-dimensional Hilbert space. 
For comparison purposes we define a uniform quantum 
distribution on the discrete version of the repeller as 



D 

-I- rep ^ 

3,3 



D 

3,3' 



(18) 



with j,/ = 1, . . . , I), and where Xi{q^p) is the finite time 
repeller defined as 



xM^p) 



if ((7,p) G Ki 
otherwise 



(19) 



with discrete time version of the repeller Ki defined as the 
set of points that remain in the allowed region for the first 
/ iterations (both for positive and negative times). The 
Irep operator can be interpreted as the coherent state 
quantization of the finite time repeller (see [l^)- Phase 
space representation for this operator, {q-,p\Irep\(l-,p)i is 
shown in top left panel of Figl6]for D = 3^ = 243. 
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Figure 5: (color online) Norm ratio (nr) for hj operators of 
the intersection quantum open baker maps with dimension 
D = 3^ = 243. Black circles, red squares and blue diamonds 
represent /c = 1, /c = 2 and /c = 5 members respectively, nr is 
also shown for corresponding Qj in dashed lines. Four Husimi 
functions {hj{q,p) = \{q,p\hj\q,p)\) for /c = 1 and /c = 5 are 
shown going form zero (blue) to maximum (red). The most 
localized cases corresponding to one periodic orbit scar (A; = 1 
and j = 2) and two periodic orbit scar (/c = 5 and j = 16) are 
shown on the right and left bottom panels respectively. Some 
cases of uniformly distributed functions on the repeller, for 
/c = 5 and j = 9 and /c = 1 and j = 18, are illustrated on left 
and right top panels respectively. Horizontal dot-dashed blue 
line is the norm ratio value obtained for the uniform quantum 
distribution on the repeller Irep- 



We have found that for the shift family the eigenstates 
corresponding to members having k > 1/2 (i.e., for which 
the openings have widths smaller or equal than V^) sud- 
denly become strongly delocalized. By strong delocal- 
ization we mean that they show nr values higher than 
the one corresponding to the Irep distribution of Eq. [181 
That is, the eigenstates are not only delocalized inside 
the discrete classical repeller but they have finite proba- 
bility in classically forbidden regions in phase space. This 
can be clearly seen in Fig. HI where the values of nr 
are shown for the shift family members k = 1,2,3, for 
D = 3^ = 243. The upper panels illustrate this remark- 
able behavior by means of phase space portraits (Husimi 
functions) of selected resonances. We should notice that 
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Figure 6: (color online) Phase space representation of 
{q,p\Irep\q,p) (left top panel) and Q32 operators for all mem- 
bers of the shift family with I = 5. Q functions are shown for 
k = 1,2,3,4,5 and j — 2^ = 32 in a D-dimensional Hilbert 
space with D = 3^ = 243. The values of the functions goes 
from zero in blue to maximum value in red. 



for D = 3^ the families at /c = 5 imply openings of "pixel" 
size and thus generate strong diffraction. The Q distri- 
bution goes asymptotically towards the uniform distribu- 
tion on the discrete version of the repeller for k < //2, 
but clearly changes its behavior above this critical value. 
In this latter case the finite probability on forbidden re- 
gions of phase space is reflected by the nr values of Q 
that saturate at about the double of that correspond- 
ing to the uniform quantum distribution on the discrete 
version of the repeller. Fig l6] shows the phase space resp- 
resentation of Irep and Qs2 for all members of shift family 
(A: = 1,2,3,4,5). 

On the other hand, the intersection family also shows 
an increase in nr values for k > 1/2 but this is now a very 
subtle effect (though a jump in nr can be seen for Q). 
All of its members have resonances distributed outside of 
the forbidden regions in phase space. They become more 
delocalized as k increases but the nr values remain always 
below the one corresponding to the uniform distribution 
on the discrete repeller. This is illustrated in Fig. [5] and 
its panels that show phase space plots (Husimi functions) 
of selected resonances. 

All this is summarized in Fig. [71 where the average 
norm ratios together with individual values for the reso- 
nances of both families are shown (in this case we have 
taken D = 3^ = 729). It is clear that a quantum phase 
transition takes place at the level of the resonances of 
the shift family. The critical value of the parameter k 
is //2, at which they become strongly delocalized and 
remain so for all values above it (we have confirm this 
result for / < 9). This applies to the long lived portion 
of the spectrum as evidenced by the two limits taken in 
the corresponding average (i.e. for the first 20 and 64 
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Figure 7: (color online) Average norm ratio < nr >, for both 
families of quantum open tri-baker maps for D = 3^ = 729. 
Top and bottom panels show shift and intersection families 
respectively. The average norm ratio for the first 20 and 64 
eigenstates are shown in red and black lines respectively, while 
circles represents individual values of norm ratios. Horizontal 
dot-dashed blue line is the norm ratio value obtained for the 
uniform quantum distribution on the repeller Irep- 



eigenstates). 

For the case of a given shift family member (fixed k > 
1/2), the localization in the forbidden area disappears 
in the semiclassical limit (increasing /), since k will be 
much smaller than the new / value. In the case of the 
intersection family there is no such phase transition but 
the delocalization of the resonances slightly grows for k > 
1/2. 

The main difference that exists for shift and intersec- 
tion families is that for one map iteration, the allowed 
regions are distributed in the whole torus for the first 
case, and only lay on the discrete time repeller for the 
second one (see Figl6]). In the quantum case, for openings 
comparable with h, the diffraction effects are very strong 
leading to notably different behaviors between both fam- 
ilies. 



IV. CONCLUSIONS 

We have investigated the behavior of the spectrum and 
the eigenstates corresponding to classical maps that share 
the same classical repeller (and other asymptotic proper- 
ties), but differ in their short time behavior. Two families 
of these maps that we call shift and intersection play a 
fundamental role in understanding their dependence on 



the shape and the area of the escape region, respectively. 
We have found that there is a quantum phase transition 
in the resonances of the shift family. For openings having 
areas smaller than h they become strongly delocalized in 
phase space, occupying forbidden regions. This makes 
difficult the quantum to classical link, even at very short 
times. Correspondingly, the spectrum oscillates in a non 
trivial way. The nature of the resonances for k > 1/2 is 
different from that corresponding to lower values of /c, as 
can be seen from Fig. [7] (upper panel). This suggests 
that the oscillation in the spectrum does not reflect an 
equivalence of the long lived sector of the extreme cases 
k = 1 and I. This will be further investigated in ^]. In 
the case of the intersection family we have found no such 
phase transition. However, there is a growing delocaliza- 
tion of the states as k grows and for k > 1/2 we have also 
identified a quick saturation of nr values for both the 
eigenstates themselves and the Q distributions towards 
the uniform distribution on the discrete repeller. In this 
family, the effect of the intersection with wider openings 
erases the probability that escapes towards forbidden re- 
gions in phase space and this prevents the phase transi- 
tion to occur here. The spectrum shrinks in this case, 
reflecting the sudden loss of probability due to a larger 
area of the opening as k grows. 

All this has very important consequences for the the- 
ory of open systems. Firstly, we mention that the fractal 
Weyl law should take into account short time dynamics 
for any given dimension D^ perhaps in its prefactor, and 
not only the invariants as the classical fractal dimension 



or escape rate. Also, the sensitivity of the spectrum with 
respect to which qubit is opened could have relevant im- 
plications to quantum computation. On the other hand, 
the strongly delocalized cases of the shift family pose an 
important question in order to understand the validity 
and performance of the semiclassical theory attempting 
to reproduce the spectrum and eigenstates of open maps. 
This is because this theory is based on highly localized 
functions, i.e. the so called scar functions which are con- 
structed with the classical information on and around 
the periodic orbits that live on the repeller. This lack 
of localization is probably due to the fact that for these 
extreme cases any quantum manifestation of the stable 
and unstable manifolds associated to the periodic orbits 
is lost. Then, it is crucial to understand how this phe- 
nomenon manifest itself even when the openings are not 
so extremely quantum in nature. This would not only 
give a better idea of how the semiclassical approximation 
works, but would tell us its validity range and perhaps 
how we can improve it. Finally, open billiards constitute 
a very interestin g fie ld of application and verification of 
our findings [15|, [l6| . Resonant micro-cavities for lasers 
can be adapted in order to realize the purely quantum 
repellers found for k > 1/2 pJ^. 

Acknowledgments 



Partial support by ANPCyT and CONICET is grate- 
fully acknowledged. 



[1] W.T. Lu, S. Sridhar and M. Zworski, Phys. Rev. Lett. 

91, 154101 (2003). 
[2] H. Schomerus and J. Tworzydlo, Phys. Rev. Lett. 93 

154102 (2004); 
[3] S. Nonnenmacher and M. Zworski, J. Phys. A 38, 10683 

(2005). 
[4] J.P. Keating, M. Novaes, S.D. Prado, and M. Sieber, 

Phys. Rev. Lett. 97 150406 (2006). 
[5] J. Wiersig, Phys. Rev. Lett. 97, 253901 (2006). 
[6] D. A. Wisniacki and G. Carlo, Phys. Rev. E 77, 

045201(R) (2008). 
[7] S. Nonnenmacher, Nonlinearity 24, R123 (2011). 
[8] J. M. Pedrosa, G. G. Carlo, D. A. Wisniacki, and L. 

Ermann, Phys. Rev. E 79, 016215 (2009). 
[9] M. Novaes, J.M. Pedrosa, D. Wisniacki, G.G. Carlo, and 

J.P Keating, Phys. Rev. E 80, 035202(R) 2009; J.M. 

Pedrosa, D. Wisniacki, G.G. Carlo, and M. Novaes, Phys. 

Rev. E 85, 036203 (2012). 
[10] L. Ermann, G.G. Carlo and M. Saraceno, Phys. Rev. 



Lett. 103, 054102 (2009). 

[11] L. Ermann, G.G. Carlo, J. Pedrosa, and M. Saraceno, 
unpublished. 

[12] L. Ermann and M. Saraceno, Phys. Rev. E 74, 046205 
(2006). 

[13] N.L. Balazs and A. Voros, Ann. Phys. 190, 1 (1989); M. 
Saraceno, Ann. Phys. 199, 37 (1990). 

[14] R.O. Vallejos and M Saraceno, J. Phys. A: Math. Gen. 
32, 7273 (1999). 

[15] C.P. Dettman, Frontiers in the Study of Chaotic Dynam- 
ical Systems with Open Problems^ (World Scientist, Lon- 
don, 2011), Chapter 11, World Scientist Series on Non- 
linear Science, Series B, Vol. 16, edited by Z. Elhadj and 
J.C. Sprott. 

[16] J. Wiersig and J. Main, Phys. Rev. E, 77, 036205 (2008); 
S. Ree and L.E. Reichl, Phys. Rev. E, 65, 055205(R) 
(2002). 



